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Abstract. Hurwitz numbers count branched covers of the Riemann sphere 
with specified ramification data, or equivalently, transitive permutation factor- 
izations in the symmetric group with specified cycle types. Monotone Hurwitz 
numbers count a restricted subset of the branched covers counted by the Hur- 
witz numbers, and have arisen in recent work on the the asymptotic expansion 
of the Harish-Chandra-Itzykson-Zuber integral. In this paper we begin a de- 
tailed study of monotone Hurwitz numbers. We prove two results that are 
reminiscent of those for classical Hurwitz numbers. The first is the mono- 
tone join-cut equation, a partial differential equation with initial conditions 
that characterizes the generating function for monotone Hurwitz numbers in 
arbitrary genus. The second is our main result, in which we give an explicit 
formula for monotone Hurwitz numbers in genus zero. 
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1. Introduction 

Hurwitz numbers count branched covers of the Riemann sphere with specified 
ramification data. They have been the subject of much mathematical interest in 
recent years, especially through the ELSV formula, given in [3], which expresses a 
Hurwitz number as a Hodge integral over the moduli space of stable curves of a 
given genus with a given number of marked points. This has led to a number of new 
proofs (see, e.g., Okounkov and Pandharipande [55] and Kazarian and Lando [TT] ) 
of Witten's conjecture [27 1 (first proved by Kontsevich [IB]), which states that a 
particular generating function for intersection numbers satisfies the KdV hierarchy 
of partial differential equations. The interest in Hurwitz numbers has much to do 
with these rich connections that their study has revealed between mathematical 
physics and algebraic geometry. There is also a connection with algebraic combi- 
natorics because of the bijection, due to Hurwitz |14j . between branched covers of 
the sphere and transitive factorizations in the symmetric group (see, e.g., Goulden, 
Jackson and Vainshtein [TT]). 

Monotone Hurwitz numbers, introduced in 9J, count a restricted subset of the 
branched covers counted by the Hurwitz numbers. The topic of [5] is the Harish- 
Chandra-Itzykson-Zuber (HCIZ) integral (see, e.g., [13], [T5], [5J]) 

I N (z,A N ,B N )= J e ^tr(A N UB N un da 

V(N) 

Here the integral is over the group of N x N complex unitary matrices against 
the normalized Haar measure, 2 is a complex parameter, and An,Bn are N x N 
complex matrices. Since U(N) is compact, In is an entire function of z G C. 
Consequently, the function 

F N {z,A N ,B N ) = N-4 ^f'^\ ^ 



is well-defined and holomorphic in a neigbourhood of z = 0, and satisfies 

I N {z,A N ,B N ) = e N2F ^ A »' B ^ 

on its domain of definition. In [S], we proved that, for two specified sequences 
of normal matrices A = (An)n=d B = C^iv)jv = i which grow in a sufficiently 
regular fashion, the derivatives F^(0,An,Bn) of Fn at z = admit an N — >• 
00 asymptotic expansion on the scale iV~ 2 whose gth coefficient is a generating 
function for the monotone double Hurwitz numbers of degree d and genus g. This 
is analogous to the well-known genus expansion of Hermitian matrix models, whose 
coefficients are generating functions enumerating graphs on surfaces (see e.g. 30]). 
In this paper, we begin a detailed study of monotone Hurwitz numbers. 

1.1. Hurwitz numbers. The single Hurwitz numbers count d-sheeted branched 
covers of the Riemann sphere by a Riemann surface where we allow arbitrary, 
but fixed, branching at one ramification point and only simple branching at other 
ramification points. Using the Hurwitz [14] encoding of a branched cover as a 
factorization in the symmetric group, we obtain the following identification with 
an enumeration question in the symmetric group: given a partition a h d and an 
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integer r > 0, the single Hurwitz number H r (a) is the number of factorizations 

(1.1) {aib 1 )(a 2 b2) ■ ■ ■ (a r b r ) = a 
in the symmetric group Sd, where 

• (a i &i), (tt2 b%), ... , (a r &r) are transpositions, 

• er is in the conjugacy class C a of permutations with cycle type a, and 

• the subgroup ((eti &i), (02 62)1 • • • , iflr b r )) ac ts transitively on the ground 
set {1,2, ...,<£}. 

Each factorization corresponds to a branched cover of the Ricmann sphere, and by 
the Riemann-Hurwitz formula, the genus g of the cover is given by the relation 

(1.2) r = d + £(a) + 2g-2, 

where 1(a) denotes the number of parts of a. Depending on the context, we will 
write H g (a) interchangeably with H r (a), using the convention that (1.2) always 
holds. 



Remark 1.1. Equation (1.1), when rewritten as (cii &i)(a2 62) • • • iflr b r )a~ = id, 
translates into a monodromy condition for the corresponding cover, in which er -1 
specifies the branching for the point with arbitrary ramification, and (ai&i), i — 
1, . . . , r specifies the (simple) branching at the remaining ramification points. The 
transitivity condition for the factorization translates to the requirement that the 
corresponding cover is connected. 

1.2. Monotone Hurwitz numbers. The monotone single Hurwitz number H r (a) 



is the number of factorizations (1.1 ) counted by the single Hurwitz number H r (a), 
but with the additional restriction that 

(1.3) 61 < 6 2 < ••• < b r , 

where <ij < 6,; by convention. As with Hurwitz numbers, depending on the context, 
we will write H g (a) interchangeably with H r (a), with the understanding that ( 1.2 ) 



holds. We will refer to a factorization (1.1) with restriction (1.3) as a monotone 
factorization of a. 

Remark 1.2. In Hurwitz's encoding, the ground set {1,2,..., d} corresponds to the 
set of sheets of the branched cover, once branch cuts have been chosen and the 
sheets have been labelled. In the case of Hurwitz numbers, the labelling of the 
sheets is immaterial, so Hurwitz numbers are usually defined to count branched 
covers with unlabelled sheets, which differs from our definition above by a factor of 
d\. However, for monotone Hurwitz numbers, the monotonicity condition depends 
on a total ordering of the sheets, so the labelling does matter in this case. Thus, for 
consistency, our convention is that both kinds of Hurwitz numbers count branched 
covers with labelled sheets. 

1.3. Main result. In |Section 4| we obtain the following theorem. This is our main 
result, and gives an explicit formula for the genus zero monotone Hurwitz numbers. 

Theorem 1.1. The genus zero monotone single Hurwitz number Hq(ol), ah d is 
given by 

1 j—i 
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where 

{2d + if = {2d + l)(2d + 2) • • • (2d + k) 
denotes a rising product with k factors, and by convention 

(2d+lf= = 

{2d + k + l)- k 

for k < 0. 

In the special case that a = (d), the partition with a single part equal to d, 
IThcorcm 1,11 becomes 

Ho((d))= =(d-l)lQ-i, 

where Cd-i — ^Cd-i) ^ s a Catalan number. This case was previously obtained by 
Gewurz and Merola [7], who used the term primitive for these fact orizations. In 



the special case that a — (l d ), the partition with all parts equal to 1, Theorem 1.1 
becomes 

H Q {{i d )) = {d-iy.2 d - 1 . 

Via the connection between monotone Hurwitz numbers and the HCIZ integral 
established in [3], this case is equivalent to a result previously obtained by Zinn- 
Justin 28J for the HCIZ integral. 

|Thcor cm l.l| is strikingly similar to the well-known explicit formula for the genus 
zero Hurwitz numbers 

(1.5) ^^(d + ^-^^n^ 



published without proof by Hurwitz [T3] in 1891 (see also Strehl and indepen- 
dently rediscovered and proved a century later by Goulden and Jackson |10j . 

There is another case in which an explicit formula similar to |Thcorcm 1.1| is 
known. This is the case where we allow arbitrary, but fixed, branching at a specified 
ramification point and arbitrary branching at all other ramification points, which 
has been studied by Bousquet-Melou and Schaeffer [5] . Given a partition ah d and 
integers r, g > 0, let Gg{a) be the number of factorizations 

(1.6) pxp2 ■ ■ ■ p r = cr 

in the symmetric group Sd which satisfy the conditions 

• Pi,P2, ■ ■ ■ ,p r € Sd, 

• a e C a , 

• (pi, P2, ■ • ■ ! Pr) acts transitively on {1, 2, ... , d}, 

and 

r 

(1.7) J2r&nk{ P] )=d + £{a)+2g-2, 

3=1 

where rank(pj) is d minus the number of cycles of pj. Each such factorization 
corresponds to a branched cover of the Riemann sphere, and (1.7), by the Ricmann- 
Hurwitz formula, specifies the genus g of the cover. 
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Remark 1.3. Note that in this case there is more freedom for the parameters a, r, 
g than for the Hurwitz and monotone Hurwitz cases above. In particular, given a 
and r, the choice for g is not unique in (1.7 1 above. This explains why we have 



used both parameters r and g in the notation Gg{a). 

Bousquet-Melou and Schaeffer [2] solved this problem in full generality for genus 
zero, using a bijective correspondence to constellations and thence to a family of 
bicoloured trees. (For more on constellations, see Lando and Zvonkin |19|.) They 
proved that 



(1.8) G » 



d\ 



|Aut a 



r((r-l)d-%) + 2) 



e(a)-2 



1(a) 

n 

3=1 



"•3 

a, 



an explicit form that is again strikingly similar to both (1.4 1 and (1.51 



The explicit formulas (1.4 1, (1.5) and (1.8) feature remarkably simple combi- 



natorial functions, but we know of no uniform bijective method to explain these 
formulas. 



1.4. Join-cut equations. The proof that we give for |Theorem 1.1[ our main result, 
involves the generating function for monotone single Hurwitz numbers 

(1.9) H(z, t, p) = J2 ^ E f E H r (a) Pa , 

d>l ' r>0 a\-d 

which is a formal power series in the indeterminates z, t and the countable set of 
indeterminates p = (pi,£>2> • • •); an d where p a denotes the product Hj=i Pay From 
a combinatorial point of view, 

• z is an exponential marker for the size d of the ground set, 

• t is an ordinary marker for the number r of transpositions, and 

• pi,p2, ■ ■ ■ are ordinary markers for the cycle lengths of a. 

In |Scction~2| we obtain the following theorem, which gives a partial differential 
equation with initial condition that uniquely specifies the generating function H. 
The proof that we give is a combinatorial join-cut analysis, and we refer to the 
partial differential equation in |Theorem 1.2| as the monotone join-cut equation. 

Theorem 1.2. The generating junction H is the unique formal power series solu- 
tion of the partial differential equation 

2t ~ <*) 2 S ^ + dlk+j + im+ > dp^dp- + ^ dp, ) 

with the initial condition [z ]H = 0. 

|Thcorcm 1.2| is again strikingly similar to the situation for the classical single 
Hurwitz numbers. To make this precise, consider the generating function for the 
classical single Hurwitz numbers 

(1.10) H ('- t .p) = E?E?E F (°)R- 

d>l ' r>0 ' ahd 



6 



I. P. GOULDEN, MATHIEU GUAY-PAQUET, AND JONATHAN NOVAK 



It is well-known (see [HUH]) that H is the unique formal power series solution of 
the partial differential equation 

(1-11) ^ = \ E («+3)PiPi 



<9H 

dp, 



iJPi+j 



d 2 H 

dpidpj 



iJPi+j 



<9H <9H 

dpt dpj 



with the initial condition [t°]H = zpi. Equation (1.11) is called the (classical) 
join-cut equation, and has exactly the same differential forms on the right-hand 



side as the monotone join-cut equation given in Theorem 1.2 There are, however 
significant differences on the left-hand side between these two versions of the join- 



cut equation. In the classical case (1.11 1, the left-hand side is a first derivative in t; 



in the monotone case (Theorem 1.2), the left-hand side is a first divided difference 



in t, and also involves differentiation in z. 

Remark 1.4. The difference in the left-hand sides between the two join-cut equa- 
tions is related to the fact that the generating functions H and H differ in the 
combinatorial role played by the indeterminate t. In the case of H, t is an expo- 
nential marker for the number r of transpositions, while in the case of H, t is an 
ordinary marker. This difference is for technical combinatorial reasons, explained 
in lScction 1.51 



Our proof of the explicit formula in |Thcorem l.l| for genus zero monotone Hurwiz 
numbers proceeds by verification using a variant of the monotone join-cut equa- 



tion given in Theorem 1.2 In general terms, this is how we obtained the explicit 
formula (1.4 1 for genus zero classical Hurwitz numbers in [10] . using the classical 
join-cut equation ( |1.11[ ). However, the technical details of this verification are quite 
different in this paper because of the change in the left-hand side between these 
two different versions of the join-cut equation. 



1.5. The group algebra of the symmetric group. In this section we express the 
generating functions H and H in terms of elements of the centre of the group algebra 
This is not an essential part of our proof of the main result, but will help 
explain why the indeterminate t for transposition factors is an exponential marker 
in H, whereas it is an ordinary marker in H. In addition, it gives a convenient 
proof of the fact that the number of monotone factorizations of a permutation a 
depends only on its conjugacy class, which we will need in |Section 2| to prove the 
monotone join-cut equation. 

First we consider the single Hurwitz numbers. For any a\- d, let C a be the formal 
sum of all elements of the conjugacy class C a , which consists of all permutations of 
cycle type a, considered as an element of Q[Sy . It is well-known that the centre of 
Q[<Sy consists precisely of linear combinations of the C Q . If we drop the transitivity 
condition for single Hurwitz numbers, the generating function for the resulting not- 
necessarily-transitive factorizations becomes 

z d if 

d>0 ' r>0 ' ahd aeC a 
z d 

= E X/ ^]E Pa I C «] C 2,l d - 2 ' 
d>0 ' r>0 ' ahd 
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where we have used the notation [A]B for the coefficient of A in the expansion of B. 
The constant term 1 corresponding to d = has been added to the generating func- 
tion t for combinatorial reasons, described as follows. When we drop the transitivity 
condition for single Hurwitz numbers, each resulting not-necessarily-transitive fac- 
torization can be split into disjoint transitive factorizations by restricting it to the 
orbits of the group ((ai b%), . . . , (a r b r )) on the ground set. Each of these orbits is a 
subset of the ground set {1, . . . , d}, and the set of transpositions that act on pairs of 
elements in a given orbit is a subset of the positions {1, . . . , r} in the factorization. 
Conversely, transitive factorizations on disjoint ground sets can be combined by 
shuffling their transpositions in any way that preserves the order of transpositions 
acting on the same component of the ground set. Thus, each factorization counted 
by t is an unordered collection of the transitive factorizations counted by H, in 
which the variables z (marking d) and t (marking r) are both exponential. From 
the Exponential Formula for exponential generating functions (see, e.g., |12|). this 
situation is captured by the equation 

H{z,t,p) = logr(z,i,p). 

Remark 1.5. The coefficient of z d t r jr\ in r(z,t,p) is in fact the image of C r 2 ld _ 2 
under the characteristic map ch of Macdonald [HJ p. 113], if one interprets the 
indeterminates Pi,P2, ■ ■ ■ as power sum symmetric functions. This can be expressed 
in the basis of Schur symmetric functions using irreducible characters of Sd, and 
then the tools of representation theory become available. While this is an interesting 
approach, we will not be using it here. 

Now we turn to monotone single Hurwitz numbers. While the monotonicity 
condition may seem artificial, it arises naturally in the group algebra Q[Sy via the 
Jucys-Murphy elements Si, defined by 

Jj = (1 i) + (2 i) H i = l,...,d. 

If we drop the transitivity condition for monotone single Hurwitz numbers, the gen- 
crating function for the resulting not-necessarily-transitive factorizations becomes 

z d 

r(M,p) = X)d7^* r &« X) W J2 lb 1 ---K 

d>0 ' r>0 ahd aeC a l<hi<-<b,.<d 
d>0 ' r>0 ahd creC a 

where h r is the rth complete symmetric polynomial. Jucys |16j showed that the 
set of symmetric polynomials in the Jucys-Murphy elements is exactly the centre 
of Q[5d], so wc obtain immediately that 

z d 

r(z,t,p) = ^7 5E *^ 5E Pq [CJMJi, • • •> Jd)- 

rf>0 ' r>0 a\-d 

This time, when we drop the transitivity condition for single monotone Hurwitz 
numbers, each resulting not-necessarily-transitive factorization can again be split 
into disjoint transitive factorizations by restricting it to the orbits of the group 
((ai bi), . . . , (a r b r )) on the ground set. Each of these orbits is a subset of the ground 
set {1, . . . , d}, and the set of transpositions that act on pairs of elements in a given 
orbit is a subset of the positions {1, . . . , r} in the factorization. However, this time, 



8 



I. P. GOULDEN, MATHIEU GUAY-PAQUET, AND JONATHAN NOVAK 



to preserve monotonicity, transitive factorizations on disjoint ground sets can be 
combined by shuffling their transpositions in only one way. Thus, each factorization 
counted by r is an unordered collection of the transitive factorizations counted by H, 
in which only the variable z (marking d) is exponential. The Exponential Formula 
for exponential generating functions then gives 

H(z,t,p) = logf(z,i,p). 

Remark 1.6. As with t(z, t, p), the coefficients of t(z, t, p) can be expressed in terms 
of the characteristic map ch of Macdonald [21 ] to provide a link with representation 
theory. This is particularly interesting in view of Okounkov and Veshik's approach 
to the representation theory of the symmetric group [23] , which features the Jucys- 
Murphy elements prominently. However, we will not be exploring this connection 
further in this paper. 

Remark 1.7. Related results on complete symmetric functions of the Jucys-Murphy 
elements have been obtained by Lassalle [2D] and Feray [B|. The recurrences they 
obtain seem to be of a completely different nature than those we obtain in this 
paper. 



1.6. Outline of paper. In Section 2 we prove the monotone join-cut equation 
of Theorem 1.2 This is based on a combinatorial join-cut analyis for monotone 
Hurwitz numbers that appears in |Section 2.1] The monotone join-cut equation itself 
is then deduced in |Scction 2.2 The join-cut analysis also yields another system of 



equations that appears in Section 2.3[ and is referred to there as a topological 
recursion. 



In Section 3 we recast the monotone join-cut equation, breaking it up into a 
separate join-cut equation for each genus, and expressing these equations in an 
algebraic form that is more convenient to solve. This is based on some algebraic 
operators that are introduced in |Scction 3.1[ and the transformed system of equa- 
tions appears in |Scction 3.2[ 



In Section 4} we prove the main result of this paper, Theorem 1.1 which gives 



an explicit formula for monotone Hurwitz numbers in genus zero. Our method is to 
repackage the formula as a generating function F, and to show that it satisfies the 
genus zero monotone join-cut equation which char acterizes the generating function 
Hq for genus zero monotone Hurwitz numbers. In 



Section 4.1 



we introduce trans- 
formed variables and use Lagrange's Implicit Function Theorem to give a closed 
form for a differential form applied to the series F. In |Section 4.2| we invert this 
differential form. In |Section 4.3[ we describe the action of the algebraic operators 
of lSection 31 on the transformed variables and deduce the main result. 



2. Join-cut analysis 

In this section, we analyze the effect of removing the last factor in a transitive 
monotone factorization, which leads to a recurrence relation for monotone single 
Hurwitz n umbers. From this, we obtain the monotone join-cut equation of |Thc-| 



1.2 



which uniquely characterizes the generating function H(z,t, p), and a 
system of equations for a different generating function that we refer to as a topo- 
logical recursion. 
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2.1. Recurrence relation. For a partition a h d, let M r (a) be the number of 
transitive monotone factorizations of a fixed but arbitrary permutation a € 5^ of 
cycle type a into r transpositions. By the discussion in |Scction 1.5| this number 
only depends on the cycle type of a, so it is well-defined, and we immediately have 

(2.1) H r (a) = \C a \M r (a). 

Theorem 2.1. The numbers M r (a) are uniquely determined by the initial condi- 
tion 

M°(a) = ! 1 ?/a = (1) ' 
I otherwise, 

and the recurrence 

(2.2) M r+1 {ai){k}) = k'm k >(a)M r (a\{k'}U{k + k'}) 

k>>l 

fe-1 



M r (aU{k',k-k'}) 
fc'=i 

fc-l r 

Y Yl MT ' ^ U { k '}) Mr ~ r ' (a\a'U{k- k'}) 



h' — l r'—O a'Ca 

for a h d, d, r > 0, k > 1. In this recurrence, mk>{ct) is the number of parts of a 
of size k' , and the last sum is over the 2^ a > subpartitions a 1 of a. 

Proof. As long as the initial condition and the recurrence relation hold, uniqueness 
follows by induction on r. The initial condition follows from the fact that for r = 
we must have cr = id, and the identity permutation is only transitive in S±. 

To show the recurrence, fix a permutation a € Sd of cycle type aU{k}, where the 
element d is in a cycle of length k, and consider a transitive monotone factorization 

(2.3) (ai b 1 )(a 2 b 2 ) ■ ■ ■ (a r b r )(a r+ i b r+1 ) = a. 

The transitivity condition forces the element d to appear in some transposition, and 
the monotonicity condition forces it to appear in the last transposition, so it must 
be that 6 r +i = d. If we move this transposition to the other side of the equation 
and set a' = a{a r+ \ b r+ \), we get the shorter monotone factorization 

(2.4) {aib 1 )(a 2 b 2 ) ■■■ (a r b r ) = a' . 



Depending on whether a r+ i is in the same cycle of a' as & r +i and whether (2.4 1 is 
still transitive, the shorter factorization falls into exactly one of the following three 
cases, corresponding to the three terms on the right-hand side of the recurrence. 
Cut: Suppose a r +\ and 6 r +i are in the same cycle of a' . Then, a is obtained 
from a' by cutting the cycle containing a r +i and 6 r +i in two parts, one 
containing a r +\ and the other containing b r+ \, so (a r +\b r+ i) is called a 
cut for a', and also for the factorization (2.3). Conversely, a r+ i and b r+ i 



are in different cycles of a, and a' is obtained from cr by joining these two 
cycles, so the transposition (a r+ i b r+ \) is called a join for a. Note that in 



the case of a cut, (2.4) is transitive if and only if (2.3) is transitive. 

For k' > 1, there are k'mk'(a) possible choices for a r +i in a cycle of a 
of length k! other than the one containing b r+1 . For each of these choices, 
(a r+ i b r+ i) is a cut and cr' has cycle type a \ {k'} U {k + k'}. Thus, the 
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number of transitive monotone factorizations of a where the last factor is 
a cut is 

k'm k >(a)M r (a\{k'}LI {k + k'}), 

k'>l 

which is the first term in the recurrence. 



Redundant join: Now suppose that (a r+1 b r+ i) is a join for a 1 and that (2.4) 



is transitive. Then, we say that (a r +i 6 r +i) is a redundant join for (2.3). 

The transposition (a r +i b r+ i) is a join for a' if and only if it is a cut 
for cr, and there are k — 1 ways of cutting the fc-cycle of a containing 6 r +i. 
Thus, the number of transitive monotone factorizations of a where the last 
factor is a redundant join is 

k-i 

M r {ali{k',k-k'}), 

k' = l 

which is the second term in the recurrence. 
Essential join: Finally, suppose that (a r +i 6 r +i) is a join for a' and that 



(2.4) is not transitive. Then, we say that (a r +i& r +i) is an essential join 



for (2.3 ). In this case, the action of the subgroup ((a\ b\), . . . , (a r b r )) must 
have exactly two orbits on the ground set, one containing a r+1 and the 
other containing b r+ \. Since transpositions acting on different orbits com- 



mute, (2.4 1 can be rearranged into a product of two transitive monotone 
factorizations on these orbits. Conversely, given a transitive monotone fac- 
torization for each orbit, this process can be reversed, and the monotonicity 
condition guarantees uniqueness of the result. 

As with redundant joins, there are k — 1 choices for a r+ i to split the 
/c-cycle of a containing Each of the other cycles of a must be in one 

of the two orbits, so there are 2^ Q ) choices for the orbit containing a r+ \. 
Thus, the number of transitive monotone factorizations of a where the last 
factor is an essential join is 

k— 1 r 

E E E M r '(a'U{k'})M r - r '(a\a'U{k-k'}), 

k' — l r' — a'Ca 

which is the third term in the recurrence. □ 

2.2. Monotone join-cut equation. Since the numbers M r (a) are a rescaled ver- 
sion of the monotone single Hurwitz numbers H r (a), we can rewrite the recurrence 



relation for M r (a) from Theorem 2.1 as a partial differential equatio n for the gen- 



Theorcm 1.2 



erating function H. The result is the monotone join-cut equation of 
which we restate here for convenience. 

Theorem 1.2. The generating function H is the unique formal power series solu- 
tion of the partial differential equation 

1 / <9H \ 1 /. , <9H .. d 2 H <9H3H\ 

2t{ Z ^- **) = 2 E (fr + J>W ^ + W ^ + ^ ) 

with the initial condition [z ]H = 0. 
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Proof. This equation can be obtained by multiplying the recurrence relation (2.2) 
by the weight 

\C a \z d+k t r p aPk 
2d\ 

and summing over all choices of d, a, k, r with d > 0, a h d, k > 1, and r > 0. The 
resulting sum can then be rewritten in terms of the generating function H via the 



defining equations (2.1 1 and (1.9 1, together with the fact that 

id- dl 

1 1 .,./'" "<"..•:<> il' 

This shows that H is indeed a solution of the partial differential equation. To see 
that the solution is unique, note that apart from d = 0, comparing the coefficient of 
z d t~ x of each side of the partial differential equation uniquely determines [z d <°]H, 
and comparing the coefficient of z d t r of each side for r > uniquely determines 
[z d t r+1 ]H in terms of [z d t r ]H. □ 

2.3. Topological recursion. In this section, we define a different type of generat- 
ing function for monotone Hurwitz numbers which is similar to the type of generat- 
ing function for classical Hurwitz numbers that has previously arisen in the physics 
literature. Specifically, by analogy with the generating function H g (xi, x%, . . . , xg) 
for Hurwitz numbers considered by Bouchard and Marino p] Equations (2.11) 
and (2.12)], consider the generating function 



(2.5) M g (x 1 ,x 2 ,...,Xi)= 



H g( a ) «! 



, Co 
ot\ ,c*2, ...,a^> 1 



1 1 j. a 2 _ 1 



where we take a — (ax, a.%, . . . , at) to be a composition, that is, an £-tuple of 
positive integers. One form of recurrence for Hurwitz numbers, expressed in terms 
of the series H g (xi,X2, ■ ■ ■ ,Xt), is referred to as topological recursion (see, e.g., [TJ 
Conjecture 2.1]; [H Remark 4.9]; 5] Definition 4.2]). The corresponding recurrence 
for monotone Hurwitz numbers, expressed in terms of the series M s (xi, X2, ■ ■ ■ , xg), 
is given in the following result. 

Theorem 2.2. For g > and I > 1, we have 

(2.6) M g (xi,x 2 , ...,xt) = SgfiS^x + XiM g -i(xi,Xi,x 2 , ■■■,xi) 

d ( x 1 M g (x 1 , ...,Xj,...x e )- XjM g (x 2 , ...,x e ) 



g 

+ ^2 ^2 x 1 M gl (x 1 ,x s )M g ^ g ,(x 1 ,Xg), 

S'=0SC{2,...,t} 

where x\, . . . , £j, . . . Xt is the list of all variables x\, . . . ,xg except Xj, xs is the 
product of all variables Xj with j G S, and S — {2, . . . ,k} \ S . 

Proof. Like the monotone join-cut equation of |Thcorcm 1.2[ this equation can be 



obtained by multiplying the recurrence (2.2) by a suitable weight and summing 



over an appropriate set of choices. In this case, the appropriate weight is 



fc-l ai-l a 2 -l . . . "f-i- 1 
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and the sum is over all positive integer choices of k, a%, a%, . . . , otg—i. In view of the 
Riemann-Hurwitz formula (1.2 1, the appropriate choice of r is 

l-i 

r = k + £ + 2g-3 + ^2aj. 

The resulting summations can then be rewritten in terms of the appropriate gen- 
erating functions by using the defining equations (2.1) and (2.5). □ 

Remark 2.1. Note that there is an asymmetry between the variable x\ and the 
variables X2, ...,X£ in (2.6), even though M. g (xi,X2, ■ ■ ■ ,x/) itself is symmetric in 
all variables. 



For small values of g and £, the recurrence (2.6) can be solved directly. In 
particular, we obtain 



Mo(ar x ) = 
M (xi,x 2 ) = 



1- Vl-4a;i 
2xi 



VI -4x^1-42:2(71 - 4^ + VI -4s 2 ) 2 ' 
If we define by y i = 1 + Xiy 2 for i > 1, then these can be rewritten as 
Mo(n) = yi, 



M (x 1 ,x 2 ) 



(yi - i)(s/2 - i)(x2 - xi) 2 



In the terminology of Eynard and Orantin [S], this seems to mean that we have 
the spectral curve y = 1 + xy 2 , but it is unclear to us what the correct notion of 
Bergmann kernel should be in our case. 



3. Intermediate forms 

In this section, we introduce some algebraic methodology that will allow us 
to solve the monotone join-cut equation. This methodology consists of a set of 
generating functions for monotone Hurwitz numbers of fixed genus, together with 
families of operators. These allow us to transform the monotone join-cut equation 
into an algebraic operator equation for these genus-specific generating functions. 

3.1. Algebraic methodology. As the first part of our algebraic methodology, we 
define three families of operators which use a new countable set of indeterminates 
x = (xi, X2, ■ ■ •), algebraically independent of p = (pi,P2, ■ ■ •)• We begin with 
lifting operators. 

Definition 3.1. Let x = (x\, x 2 , ■ ■ ■) and p = (pi,P2,---) be countable sets of 
indeterminates. The ith lifting operator Aj is the Q[[x]] -linear differential operator 
on the ring Q[[x, p]] defined by 

d 

Ai = Vfc^ — , i>l. 
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The combinatorial effect of Aj , when applied to a generating function, is to pick 
a cycle marked by p^ in all possible ways and mark it by kx\ instead, that is, by 
x\ once for each element of the cycle. Note that AiXj = for all j, so that 



A? = £ jk 

j,k>i 



j + k 
X J - 





d 2 

dpjdpk ' 



Accompanying these lifting operators, we also have projection operators. 

Definition 3.2. Let x = (xi,x 2 ,--) and p = (pi,P2,---) be countable sets of 
indeterminates. The zth projection operator IT is the Q[[p]]-linear idempotent 
operator on the ring Q[[x, p]] defined by 

n, = [^] + 5> fc k fe ], »>i. 

k>l 

The combinatorial effect of IT, when applied to a generating function, is to take 
any cycle marked by x\ and mark it by pk instead. The combined effect of a lift 
and a projection when applied to a generating function in Q[[p]] is given by 

Finally, we introduce splitting operators. 

Definition 3.3. Let F(xi) be an element of Q[[x, p]], considered as a power series 
in Xi, and let j > 1 be an index other than i > 1. Then the i-to-j splitting operator 
is defined by 

Split F( Xi ) = ^ilM^llM. +F ( ), 



so that 



Splits? - "'' " 2 - 2 ■■ J : 



X^ Xj ~\~ X^ Xj ~ p ' ■ ' — p X^Xj 



Combinatorially, the effect of Split^j on a generating function is to take the 
cycle marked by X{ and split it in two cycles, marked by Xi and Xj respectively, 
in all possible ways. The combined effect of a lift, a split and a projection on a 
generating function in Q[[p]] is 

ILITa Split Ai = V {i+j)piPj^ . 

n-2 dp t+j 

As the second part of our algebraic methodology, we define the generating func- 
tions 

(3-D fi 9 = ££%^> 9>0. 

d>la\-d 

Thus, Hg is the generating function for genus g monotone single Hurwitz numbers 
where, combinatorially, P\,Pii ■ ■ ■ are ordinary markers for the parts of a, and there 
is an implicit exponential marker for the size d of the ground set. 
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3.2. Recasting the monotone join-cut equation. We are now able to recast 
monotone join-cut equation as an algebraic operator equation involving our genus- 
specific generating series H g and the three operators that we have introduced above. 

Theorem 3.4. 

(i) The generating function A{Rq is the unique formal power series solution of 
the partial differential equation 

(3.2) 



AxHo = n 2 Split AiH + (AiH ) + n 

l->2 



with the initial condition [po^iJAiHo = 0. 
(ii) For g > I, the generating function AiH g is uniquely determined in terms of 
AiHq, AiHi, . . . , AiH g _i by the equation 



(3.3) 



1 - 2A!H - n 2 Split ) AiHg = A^Hg.! 



2 / 



9-1 

E 

9' = 1 



AiH„' AiH„ 



(Hi) For g > 0, the generating function H g is uniquely determined by the generating 
function AiH s and the fact that [po]H g = 0. 

Proof. As with |Theorem O and Theorem 2.2| this result is obtained by multiplying 
the recurrence from Theorem 2.f by a suitable weight and summing over a set of 
possible choices. Consider the generating function 



g>0 g>0 d>la\-d 



Hg(a)p c 



where u is an ordinary for the genus g. In view of ( |1.2[ ) and (2.1), we have 

M d +^+ 2 9- 2 (a)p a 



F 



g>Q 



g>0 



'EEE 

d>0 a\-d k>l 



M d+k+l( a )+2g-l( a u | fc |) paX k 



n 3 >if (a S(«)i 



Thus, by multiplying the recurrence (2.2) by the weight 



u 9 p a x\ 



Ilj>iJ m,(a Wa)! 

and summing over all choices of g,d,a,k with g > 0, d > 0, a \~ d, k > 1, and 
r = d + k + £(a) + 2g — 2, we obtain the partial differential equation 

(3.4) A X F -Xi=IL 2 Split A X F + uAfF + (A X F) 2 . 

1^2 

To show that A\F is the unique solution of this partial differential equation with 



[u°pqXi]AiF = 0, note that (3.4) exactly captures the recurrence of Theorem 2.1 
so each non-constant coefficient of AiF is uniquely determined. 

Extracting the coefficient of u 9 from ( 3.4 ) and rearranging terms gives the mono- 
tone join-cut equations of the theorem statement for 5 = and g > 1. 

Finally, note that given AiH g , we can compute 



n 1 A 1 n q 



VV 



H g (a)p a 



(d-l)\ 

d>lahd y ' 
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which uniquely determines every coefficient of H ff except for the constant term. □ 

Remark 3.1. This form of the monotone join-cut equation is technically slightly 
stronger than the one given in |Theorem 1.2\ since it is obtained from the recurrence 
relation in Theorem 2.1 by using a less symmetric weight. 

Remark 3.2. The monotone join-cut equation for higher genera will be the subject 
of a forthcoming paper 18 , which we won't discuss further here. 

4. Transformed variables and proof of the main result 



In this section, we prove the main result, Theorem 1.1 Our strategy is to define 
the series F by 

d>l ahd 1 1 3=1 

and then to show that the series AiF satisfies the genus zero monotone join-cut 
equation (3.2 I. 

Remark 4.1. We initially conjectured this formula for genus zero monotone Hur- 
witz numbers after generating extensive numerical data, using the group algebra 
approach described in Section 1.5 together with the character theory and gener- 
ating series capabilities of Sage |25j . In particular, the case where a has t(ot) = 3 
parts was very suggestive, since the formula then breaks down into a product of 
three terms. This was also our first indication of the striking similarities between 
monotone Hurwitz numbers and classical Hurwitz numbers. 

4.1. Transformed variables and Lagrange inversion. In working with the 
series F, it is convenient to change variables from p = (pi,P2, . . .) to q = (qi, q2, ■ . .), 
where 

-2j 

'2k\ 



(4-2) q = Pj I I ' i — 1- 

This change of variables is invertible, and can be carried out using the Lagrange 
Implicit Function Theorem in many variables (see 12J). The first result expressing 
F in terms of the new indeterminates q involves the differential operator 

(4.3) P = E^ 

fc>i 1 K 

Theorem 4.1. Let 7 = J2k>i Ck)l k ant ^ = 

E fc > 1 (2fc + l)( 2 fc fe )<z fc . Then 
(1 -7) 3 

(2V - 2)<2V - 1)(2£>)F = ^ — - 1. 

1 - 77 

Proof. From (4.1 ), for any ahd with d > 1, we have 

[Pa] (2V - 2)(2V - 1)(2P )F = ^ (2d - 2)^) l| 

{-l)H{a)\ (2-2d\ ^ (2 aj 



|Auta| V 1(a) 



n 

3=1 
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and we conclude that 

(4.4) [p Q ](2Z> - 2)(22? - 1)(22?)F = [q a ](l - jf~ 2d . 



Now let <f>j = (1 — 7)" J , so that (4.2) becomes qj = Pj<fij, j > 1. Then, from 
the multivariate Lagrange Implicit Function Theorem [T5J Theorem 1.2.9], for any 
formal power series <E> € Q[[q]], we obtain 

d 



[Pa]$ = [&*]$&» det £ 



% ^— log 4 



= [<?<*]$ </> a det ^ - - 
Then we have 



7 V J 

\-2d 



»,3>1 



where Q = rii'>i Then we have (f> a = (1 — 7) 2d , and using the fact that 
det(7 + M) = 1 + tr(M) for any matrix M of rank zero or one, we can evaluate the 
determinant as 



d 

<1H I Sij - qj — logi 



= 1 - 



2kq k 

^ I-7 

fe>i ' 



1 — r/ 
1-7' 



Substituting, we obtain 



[Pa]^ = [9a] 



(1 - 7?)$ 



1 (1 _ 7 )2d+l ' 

Comparing this result with (|4.4[), we obtain 



ba](2Z> - 2) (22? - 1)(22?)F = [ Pa ] 



(l- 7 ) ; 



1 



V 



for a h d and d > 1, and computing the constant term separately, the result follows 
immediately. □ 

4.2. Inverting differential operators in the transformed variables. In order 
to use Theorem 4.1 to evaluate AiF, we need to invert the differential operators 
TD — 2, TD — 1, and 22?. We will work with the transformed variables q, and thus 
introduce the additional differential operators 



D» 



Pk 



dp k 



Ei 



Qk 



dq k 



^kq k 

k>l 



d 



As Q-linear operators, the operators Di,D 2 ,... and 2? have an eigenbasis given 
by the set {p a - ot h d, d > 0}, and consequently they commute with each other. 
Similarly the operators E!,E 2 , . . . and £ have the set {q a : a h d, d > 0} as an 
eigenbasis and commute with each other. However, these two families of operators 
don't commute with each other. By using the relation (4.2) to compute the action 
of Efc on pj , we can verify the operator identity 



Ei. = Di 



2q k 



1 



It follows that 



1-71 



V. 



2?, 



k > 1. 



and we can deduce the identity 
(4.5) D fe = E fe 



2gfc 

1 — T] 



£. 



k > 1. 
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Thus, we can express these differential operators for p and q in terms of each 
other. In the following result, we apply these expressions to invert the differential 
operators that appear in |Theorem 4.1| 

Theorem 4.2. For k > 1, we have 

DfcF = 2fc(2ibl) ( *> - g 2{3 + k)(2k-l) (j) 
Proof. As notation local to this proof, let 



F'" = (22? -2) (22?- 1)(2D)F, F" = (22? - 1)(22?)F, F' = (22?)F. 
To prove the result, we use the operator identity 



(4.6) 



(i - 7 r(2£ - o ((i - 7 r G ) = r Ji(2D- o(g), 



which holds for any integer i and any formal power series G. This allows us to 
express the differential operators 22? — 2, 22? — 1, and 22? in terms of the operators 
25 — 2, 25 — 1, and 2£ 1 which we can invert by recalling that they have {q a : a h 
d, d > 0} as an eigenbasis. 

We proceed in a number of stages. First we invert 22? — 2 by applying (4.6 ) with 
i = 2 to Theorem 4.1 obtaining 

F" = (22?-2)- 1 (F"') 
= i + (22?-2)- 1 



1 — Tj 



i + (l- 7 ) 2 (2£- 2)" 1 (1) 
i i / 



= £-i(i-7)', 

after checking separately that [pi]F" = 2. (We need to check this since the kernel 
of 22? — 2 is spanned by p\.) 

Next we apply to F" via (4.5 1. This is straightforward, and gives 



D fe F" = 



(1 - 7 ) 2 (2k\ 



1 — 77 



qk- 



Now we invert 22? — 1 by applying (4.6 1 with i = 1, which gives 



D fc F' = (22?- I)- 1 (D fc F" 
= (l- 7 )(2£-l)- 1 
1 



2fc\ 
k 



= (l- 7 ) 



2k -1 



2k\ 
k 



qk- 
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Finally, we invert 21? by applying (4.6) with i = 0, giving 
D fc F = (2V)- 1 (B k F') 

= < 2 ^ , (< 1 -"»s L t(T)" 




2j + 1 (2j 



2k(2k - 1 



2k-l\j 
^2(j + fc)(2fc-l)Vj 




Again, the constant term needs to be checked separately, since the kernel of (22?) 
consists of the constants, but clearly D^-F has no constant term. □ 

4.3. The generating function for genus zero. In order to work consistently in 
the tranformed variables q, it will be useful to have descriptions of the projection 
and splitting operators in terms of q. When considering these operators, the change 
of variables from p to q also corresponds to a change of variables from x to a new 
countable set of indeterminates y = (yi, j/2, • • where we impose the relations 



(4.7) 



Ui = Xi(l - 7)" 



We can express the indeterminates p and q in terms of each other using (4.2), so 
we can identify the rings Q[[p]] and Q[[q]]. Since (1 — j)~ 2 is an invcrtiblc clement 
in this ring, we can further identify the rings Q[[p,x]] and Q[[q, y]] using (4.7). In 
this bigger ring, we have the operator identities 



n, 



Split G( Xi ) = 



E 

fc>i 



Pk[Xi 



fe>i 



Xj (^Xi ^ X i G {jX j j 



. G(0) = »g(«i) - wcfo) +g(0) 

Vi - Vj 



so the projection and splitting operators are just as easy to use with either set of 
indeterminates. 

Remark 4.2. For completeness, note that the lifting operators can also be described 
in terms of q and y, although the expressions are somewhat more complicated. That 
is, using (4.5), we obtain the expression 

k>l LK k>l 



A; 



4^(1-4^)- 

(1-7?) 



k>l 



d_ 
dq k 



dpk 'dqk, 
We are now able to evaluate AiF in the indeterminates q and y. 
Corollary 4.3. We have 



Vk 



dyk 
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Proof. From [Theorem 4.2[ we have 



AxF = V ^D fc F 



E 



i 



2k 



(2j + l)k [2j\[2k\ k 
1 U Vfc 



k - 2(2k - 1) V k J" 1 2(j + k)(2k - 1) V j J\ k 

= n 2 (2G(y 1 ,0)-G(y 1) y 2 )), 

where the power series G(yi,y 2 ) is defined by 



j>0 k>l 



Then, the computation 



G(yi,y 2 )= / yit(l-4 yi t)-?(l-4y 2 t) 



2 

t 



-Hi 



2(2/1 - 2/ 2 ) 

V -± 

2(2/i - 2/2) 



(l-4»it)*(l-4ifct)-3 
l-42/i 



1 - 



1-42/2 



completes the proof. 



□ 



In the following result, using the above explicit expression for A X F, we uniquely 
identify F as the generating function for monotone single Hurwitz numbers in genus 
zero. 

Theorem 4.4. The series F satisfies the genus zero monotone join-cut equation 

AiF = n 2 Split AiF + (A1F) 2 + x x . 

1^2 

Thus, F = H is the generating function for monotone single Hurwitz numbers in 
genus zero. 



Proof. Corollary 4.3 gives AiF = U 2 A(yi, y 2 ), where 



-4(2/1,2/2) = 1-^1-42/1 



We need to check that the expression 



2/1 



2(2/1 - 2/2) 



1-42/! 
1-42/2 



AiF - n 2 Split AiF - (AiF) - x x 

l-s-2 
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is zero. To do so, we rewrite each of the terms in this expression as 

'2/ 2^(2/1,2/3) - 2/1^(2/2,2/3) 
2/i - 2/2 



AiF 


= n 2 n 3 


n 2 Split A : F 


= n 2 n 3 


i-s-2 




(A^) 2 


= n 2 n 3 


xi = 2/i(l - 7 ) 2 


= n 2 n 3 



n 2 n 3 \A(yx,y 2 )A(y u y 3 



1 W_ 1 



to get an expression of the form 

n 2 n 3 5 (2/1,2/2,2/3) ■ 

The series B(y±, 2/ 2 , 2/3) itself is not zero, but a straightforward computation shows 
that the series 

§£(2/1,2/2,2/3) + ^ B (yi>y3,y2), 

obtained by symmetrizing with respect to 2/2 and 2/3, is zero. Thus we have 

n 2 n 3 B{ yi , 2/ 2 , 2/3) = n 2 n 3 (^(2/1, 2/2, 2/3) + §£(2/1, 2/3, 2/2)) = 0, 

which completes the verification. The fact that F = H follows immediately from 
ITheorem 3.4l □ 

Finally, we are now able to deduce our main result. 



Proof of \Theorem In view of (3.1 1 with g = 0, applying Theorem 4.4 gives 



H (a) =dl[ Pa ]F 



for any partition a of d > 1. The result follows immediately from (4.1 1. □ 
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